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ABSTRACT 

We report the identification of a new cataclysmic variable (denoted as CV2) and a probable 
I ■ microlensing event in the field of the globular cluster M22. Two outbursts were observed for 

Q ' CV2. During one of them superhumps with =0.08875 d were present in the light curve. 

?H ' CV2 has an X-ray counterpart detected by XMM-Newton. A very likely microlensing event 

jy^ i at a radius of 2'3 from the cluster center was detected. It had an amplitude of AV = 0.75 mag 

j^j , and a characteristic time of 15.9 days. Based on model considerations we show that the most 

likely configuration has the source in the Galactic bulge with the lens in the cluster. Two 
outbursts were observed for the already known dwarf nova CV1. 

, Stars: dwarf novae - novae, cataclysmic variables - globular clusters: individual: (M22, 



5^ ■ NGC 6656) - Gravitational lensing 

1 Introduction 



The dense central regions of globular clusters (GCs) offer a unique opportunity 
for studies of dynamical processes in stellar systems. It is well established that 
even a small population of close binaries can drive the evolution of the entire 
cluster (Hut et al. 1992, 2003). It is expected on theoretical grounds that large 
numbers of cataclysmic variables (CVs) should be present in the central parts 
of GCs (Fabian et al. 1975; Hut and Paczyhski 1984). Several tens of candidate 
cluster CVs have been reported over the last few years based on observations col- 
lected with the Chandra (e.g., Hannikainen et al. 2005, Heinke et al. 2005) and 
XMM-Newton (e.g., Gendre et al. 2003; Webb et al. 2004) telescopes. However, 
only a few CVs in GCs have been confirmed with spectroscopic observations. 
Similarly, there are only a few objects for which dwarf nova (DN) type outbursts 
have been observed (see Kaluzny et al. 2005 for more extended summary of the 
subject). 

This contribution continues a series of papers devoted to the search for CVs 
in GCs based on the rich photometric data base collected by the CASE collab- 
oration (Kaluzny et al. 2005). We report here the results for the cluster M22 
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Table 1: Positions of eruptive objects observed in M22 



Name 


RA(2000.0) 


Dec(2000.0) 


Distance 








from center 


M22 CV1 


18 h 36 m 24?66 


-23°54'35"5 


0!40 


M22 CV2 


18 h 36 m 02?72 


-23°55'24"6 


5'.51 


M22 microlens 


18 h 36 m 22?40 


-23°56'29"4 


2!33 



(NGC 6656) which is located about one-third of the way between the Sun and 
the Galactic bulge. Two DN type outbursts were detected for the previously 
known variable CV1 and for the newly identified variable CV2. We have also 
observed a very likely microlensing event caused by a cluster star. 

2 Observations and Search for Erupting Objects 

The CASE project is conducted at Las Campanas Observatory. For the survey 
we used the 1.0-m Swope telescope equipped with a 2048 x 3150 pixel SITE3 
CCD camera. With a scale of 0.435 arcsec/pixel the field of view was 14.8 x 
23 arcmin 2 . A fraction of the images of M22 was taken in a subrastered mode 
with a field of 14.8 x 15.6 arcmin 2 . The analysis presented in this paper is based 
on a uniform set of trimmed images with a field of view 14.8 x 11.6 arcmin 2 
(longer axis along E-W direction). The cluster core was located approximately 
at the center of the images. 

The cluster was monitored during the 2000 and 2001 seasons. A total of 
2006 images were taken through the V filter with exposure times 90-300s. In 
addition, we have obtained 390 images in the B filter with exposure times 140- 
300s. The number of V-band frames taken per night ranged from a few to 62. 
The median seeing was 1"41 and 1"43 for the V and B bands, respectively. 

The photometry was extracted with the help of the Difference Image Anal- 
ysis Package (DIAPL) written by Wozniak (2000) and recently modified by W. 
Pych. The package is an implementation of the method developed by Alard and 
Lupton (1998). The procedure used to search for possible " out bursting" objects 
is described in detail in Kaluzny et al. (2005). The search was based on the V 
band data collected on 71 nights for which at least one image with seeing better 
than 1"6 was available. A reference frame for the V filter was constructed by 
combining 19 individual frames taken on the night of 2000 Sep 2/3. The seeing 
for that frame was FWHM = 1"03. Profile photometry, as well as aperture cor- 
rections for the reference images, was extracted with DAOPHOT/ALLSTAR 
(Stetson 1987). These measurements were subsequently used to transform the 
light curves from differential flux units into instrumental magnitudes. The final 
transformations onto the standard BV system were based on observations of 
several Landolt (1992) fields which were collected on one of the nights during 
which M22 was observed. 

Searches for erupting objects in M22 led to the detection of three objects. 
In addition to the previously known DN CV1 (Anderson et al. 2003; Bond et 
al. 2005), we identified another DN (which we provisionally name CV2) and a 
probable microlensing event. The equatorial coordinates of these three objects, 
as well as their angular distances from the cluster center (a 2 ooo = 18 /l 36 m 24?2, 
£2000 = -23° 54' 12", Harris 1996), are listed in Table 1. 
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3 Outbursts of CV1 

The variability of CV1 was first noted by Sahu et al. (2001). They interpreted 
the May 1999 brightening episode detected on HST/WFPC2 images as a mi- 
crolensing event. Subsequently Anderson et al. (2003) argued that in fact the 
variable is a CV and that the observed event was due to a DN type outburst. 
The outburst lasted about 25 days with an approximate amplitude of 3 mag, 
peaking at J« 15 mag. Two other outburst episodes, one in 2002 and one in 
2003, were reported by Bond et al. (2005). For the event in 2002 only one expo- 
sure was obtained showing the variable at 7=15.1. During the 2003 episode the 
variable was observed at I m ax ~ 15.7 mag, but the actual peak of the outburst 
was very likely not observed. 
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Fig. 1. Light curves of CV1: full span in the 2000-2001 seasons (upper panel), during the 
2000 Aug/Sep outburst (middle panel) and nightly segmented during the outbursts (lower 
panels). In the segments small ticks are every 0.05 days. 



As can be seen in Fig. 1 another two outbursts of CV1 were detected in our 
data. The first one occurred on August/September 2000 and lasted approxi- 
mately 20 days. The observed maximum reached V = 16.75 mag, but we failed 
to cover the actual peak of the outburst. There is no evidence for the presence 
of superhumps in the light curve. The second outburst episode of CV1 in our 
data set occurred in 2001. On June 3, 2001 the variable had a mean magnitude 
of V = 17.22. Unfortunately, we observed the cluster on only one night during 
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this event, and observations collected seven days later showed that CV1 had 
already faded to a low state. 

Our data taken together with the data of Bond et al. (2005) cover the observ- 
ing seasons 2000-2004. A total of 4 outbursts were observed during that period 
suggesting that the average recurrence time is rather long, probably exceeding 
150 days. The duration of the three outbursts with reasonably complete time 
coverage ranges from about 20 to about 25 days (Sahu et al. 2001; Bond 2005; 
this paper). 

In Fig. 2 the location of CV1 during the 2000 eruption is marked on a 
color-magnitude diagram of M22. The cluster main sequence shows a noticeable 
spread due to the presence of significant differential reddening (0.34<i?(i?— V)< 
0.42 (Richtcr et al. 1999)). In the following analysis we use £(B-T/)=0.38±0.04 
for the cluster reddening. The variable had a mean color of <B— F>=0.46±0.09 
near maximum light which implies an unreddened color of (B — V)q = 0.08± 
0.13. This is within the range of colors exhibited by DNe during outbursts 
(Warner 1995). We were unable to obtain reliable photometry of the variable 
on many images collected in the low state due to crowding effects. From archival 
HST/WFPC2 images (GO 8174 program) observed on 2000 June 28 we obtained 
V m i n = 18.96. The photometry was extracted with the help of HSTphot package 
(Dolphin 2000a, b). In Fig. 1 we adopt this value of V m %n to mark all out-of- 
outburst observations of CV1. 
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Fig. 2. Locations of observed cataclysmic variables and the microlensing event on the color- 
magnitude diagram of M22 
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4 Variable CV2 

Fig. 3 shows the light curve of the newly detected candidate DN. A fading 
branch of an apparent outburst was observed during the period 2000 July 28- 
August 3. During that time the brightness declined from V~17.7 mag to V~19.0 
mag. The main outburst was followed by an "echo" eruption lasting about two 
days. The light curve for the first three nights of observations covering the out- 
burst (lower panel in Fig. 3) shows the presence of supcrhumps. Their shape and 
full amplitudes of «0.18 mag are typical for SU UMa type DN (Warner 1985). 
After subtracting nightly means from the light curve we calculated the power 
spectrum with ANOVA statistics (Schwarzenberg-Czerny 1996). The first two 
harmonics of the Fourier series were used in the fit. The resulting periodogram 
is shown in the upper panel of Fig. 4. The highest peak corresponds to a su- 
perhump period of P sh = 0.08875 ±0.00012 days (127.80 ±0.18 min), while the 
second peak is its subharmonic. The lower panel of Fig. 4 presents the light 
curve folded with P s h and averaged in 0.05 phase bins. The orbital period of 
CV2 can be estimated using the empirical relation between P s h and P or b first 
noted by Stolz and Schoembs (1984). Using the calibration derived by Skill- 
man and Patterson (1993) we obtain P or b~ 0.08496 d (122.3 min). This places 
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Fig. 3. Light curves of CV2: full span in the 2000-2001 seasons (upper panel), during the 
2000 July/August eruption (middle panel) and nightly segmented during the superoutburst 
(lower panels). Note the presence of an "echo" outburst just after the superoutburst. 
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Fig. 4. Upper panel: the ANOVA power spectrum of CV2 of the JD=2451754-56 light curve 
of CV2. Lower panel: the light curve folded at frequency 11.267 cycles/day. 



the variable at the lower edge of the gap observed in the distribution of orbital 
periods for CVs (Ritter and Kolb 2003). 

Figure 5 presents finding charts for CV2 showing the variable in quiescence 
and in outburst. The variable is plotted on the cluster color- magnitude diagram 
in Fig. 2. We derive B — V = 0.95 ± 0.05 mag and B — V = 0.45 ± 0.05 mag for 
the low state and at maximum light, respectively. Assuming a reddening of 
E{B -V) = 0.38 ± 0.04 (Richter et al. 1999) we estimate the unreddened color 
at maximum light to be (B — V)o s=s 0.07 mag. Most DNe at maximum light 
have unreddened colors in the range (B — V)o = 0.0± 0.10 mag (Warner 1976). 
The variable was observed at V « 19.0 mag in quiescence. This brightness is 
consistent with the hypothesis that CV2 belongs to M22, despite being located 
at 3.9 core radii from the cluster center. Assuming cluster membership the 
variable has My«5.5 mag in quiescence (we adopt (to— M)y=13.5 mag for M22, 
Harris 1996). This places CV2 in the range of absolute magnitudes observed for 
bright non-magnetic DN. 

Further confirmation of the cataclysmic nature of CV2 comes from the X- 
ray observations of M22 obtained with XMM-Newton. Webb et al. (2004) pub- 
lished a list of 50 X-ray sources detected in the 30' x 30' field centered on the 
cluster. One of the sources, object #40, is a very probable counterpart of 
CV2. The equatorial coordinates of the X-ray object are a^ooo = 18 /l 36 m 02?96, 
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Fig. 5. Finding charts for CV2 showing the variable in quiescence (left) and during outburst 
(right). North is up and East is to the left. The field of view is 13" on a side. 

^2000 = — 23 55 26"42, offset from the optical coordinates of CV2 listed in Ta- 
ble 1 by (Aa = 3"6, A6 = 1"8), well within the 5" uncertainty circle of the 
X-ray source position. The XMM-Newton observations were obtained on Sep 
19-20, 2000. That was 20 days after our last observation in the 2000 season and 
45 days after the eruption of CV2. The measured count rate in the 0.5-10.0 keV 
band for source #40 was 2.74 times lower than that measured for source #36 
associated with CV1. For CV1 Webb et al. (2004) give an unabsorbed flux of 
9.8 x 10~ 14 erg cm _2 s _1 . Assuming that both CVs belong to the cluster at a 
distance of 3.2 kpc (Monaco et al. 2004) we can estimate the X-ray luminosity 
of CV2 during quiescence as Lx ~4.4 x 10 31 erg s _1 . This value is consistent 
with X-ray luminosities observed for field non-magnetic CVs (e.g., Baskill et al. 
2005). 

5 Microlensing Event in M22 

Our search for erupting objects in M22 also led to the detection of a source 
which underwent one episode of increased luminosity. Its light curve is shown in 
Fig. 6. The source is located 140 from the cluster center in a relatively crowded 
region (see the chart in Fig. 7). The brightness of the object increased by w0.7 
mag over 20 days. Around 2000 August 5 it reached a maximum of VPS19.1 and 
then faded to a constant level of V ~ 19.9 mag. On the cluster color-magnitude 
diagram (Fig. 2) the star is placed among the cluster main sequence stars. 
However, we note that at minimum brightness it is only marginally detectable 
on B images leading to a large uncertainty in the measured color. The variable 
is relatively red at the at the maximum light with B~V = 0.97 ±0.09 mag and 
(B — V)o~0.59 mag. * This excluded possibility that the object belongs to DNs 
(Warner 1976). At minimum light we measured B — V = 1.24 ± 0.22 mag what 
is consistent with a lack of any color change during an outburst. 

The observed behaviour of the star is characteristic of a gravitational lens 
event. To check this hypothesis we use a 5-parameter model of a single lens 
event and fit a light curve of the standard form (Paczynski 1986): 

F(t) = A(u(t)) F s + F h 

*According to Schlegel et al. (1998) the total reddening in the cluster direction amounts 
to E(B — V) = 0.33 mag. That implies that any stars located in the cluster field cannot be 
significantly more reddened than the cluster itself. 
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Fig. 6. Full span (top) and event (bottom) light curves of the probable microlcns in M22. 




MVM +4 

where F s is the (unmagnificd) flux of the star undergoing lcnsing, F\> - the 
blended flux, to - the epoch of maximum, - a characteristic (Einstein) time 
of the event, and u the impact parameter expressed in units of the Einstein 
radius. The Einstein radius for a lens of mass M located at a distance cJl, and 
a source at a larger distance ds is: 



4GMd L (d s -d L ) 

rE= i~ — d S — (2) 

where G is the gravity constant and c the speed of light. 
We transform the measured stellar magnitudes Vi and their error estimates 
AVi into flux units (Fi, ai) and fit the model minimizing x 2 '- 

X 2 =E ( ^ + f b ~ F ' )2 (3) 



where Ai = A(u(ti)). 
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Fig. 7. Finding chart for the probable microlens in M22. North is up and East is to the left. 
The field of view is 13 " on a side. 

The fitted parameters, with source and blend fluxes converted back into 
stellar magnitudes V S} Vb are: 

to = 2451759.70^11 t E = 15. d 9±^ Uo = 0.541^ 

V B = 19.92tg;{g Vb = 24.8l2° 

X 2 = 127 DOF=1923 

The fit, shown in Fig. 8 is not well constrained due to the Wozniak and 
Paczyhski (1997) degeneracy. The fitted parameters are highly correlated within 
the confidence regions in parameter space. The event can be fitted by a relatively 
shorter lasting and weaker amplification of a brighter source with a fainter blend 
as well as by a longer lasting and stronger amplification of a fainter source with 
a brighter blend. 

The source and the lens can be located in the bulge, in M22, or anywhere 
along the line of sight. Since the large majority of objects of interest belong 
either to M22 or to the bulge, we shall consider primarily these locations. First 
we assume that the lens belongs to the cluster and compare the a priori relative 
probability of the source being in the bulge or in the cluster. 

We adopt the distances cZm22 = 3.2kpc to the cluster (Monaco et al. 2004) 
and c?b = 8.0kpc to the bulge. The Einstein ring projected into the source plane 
has a radius r^ = r-E dg/d^. Only sources within from the observer - lens line 
are viable candidates for lensing. 

For a rough estimation of the relative lensing probability we adopt a sim- 
plified model of the cluster assuming spherical symmetry and a density profile 
n(r) ~ r~ 3 outside the core. The line of sight to the lensed source passes at 
&M22 = 2.17pc from the cluster center while the core radius r c = 1.30pc (Trager 
et al. 1995). In our model the surface density of cluster stars is related to its 
space density along the line of sight: 

r+oo b 3 

Nm22 = J-co n >2 + x 2y3/2 dx = 2n b (4) 

where no is the 3D star density at a distance b from the cluster center. The 
probability of finding a cluster source within the Einstein radius from a lens of 
mass M is: 



M22 : 



1 f +0 ° f + °° f 4GM 

dx n(x) / dy< n(y) it — = — x 



^M22 



(d M22 + y)(y-x) ^ NM J_^ bM22 (5) 

CIM22 +X I C 2 
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Fig. 8. The standard (5 parameter) binary lens model of the microlensing event. The data 
points were nightly binned for clarity of the plot. 



where x and y are the lens and source positions respectively, measured along 
the line of sight. We approximate (c?M22 + 2/)/(^M22 + x) w 1 and then solve the 
double integral analytically. 

For a source in the bulge the lens position within the cluster is irrelevant, 
and one gets a probability of lensing: 



AttGM d B (d B -dM22) 



(6) 



- d M 22 

The relative probability of the lensed source belonging to the cluster is 
flvr22/-PB=9.04x 10 -5 7Vm22/-/Vb~5x 10~ 4 where we use estimates of the surface 
densities of stars belonging to M22 and the bulge from Albrow et al. (2002). 

Similarly, assuming that the source is located in the bulge, one can check 
the relative probability of finding the lens in the bulge or in the cluster. The 
calculations are analogous to those presented above, but the surface density of 
sources should be replaced by surface mass density of lenses. The line of sight 
passes at w 12° or & B ~ l-7kpc from the bulge center. Calculations give 



(7) 
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47tGSm22 d,M22{d>B — ^M22) ,„s 
C 2 d B 

where Sb, £m22 stand for the surface mass density in the bulge and in M22, 
respectively, measured along the line of sight. The geometrical factors are of the 
same order, so Pb/Pm22 = 0.42£b/£m22- Assuming that surface mass densities 
are proportional to the surface densities of stars we get Pb/Pm22 ~0.09. The 
models of mass distribution in M22 and bulge suggest even an smaller value. 
Thus our estimates show that the location of the source in the bulge with lens 
in M22 is the most likely situation. The location of both objects in the bulge is 
> 11 times less likely, while the location of both objects in M22 is w 2000 times 
less likely. In the following we assume that the source is in the bulge and the 
lens in the cluster. 

The proper motion of the lens relative to the source is dominated by the 
motion of the cluster as a whole (Gaudi 2002), yU ro i = 10.9 mas y _1 (Peterson 
and Cudworth 1994). The velocity dispersion of stars in M22, a = 11.4 km s _1 
(ibid.) corresponds to fi = 0.75 mas y _1 , and the motion of stars in the bulge 
with velocity dispersion components (oi,<7b) = (93,79) km s _1 (Han and Gould 
1996) corresponds to /i?s2 mas y _1 , so they can be neglected to simplify the 
analysis. By the definition of Einstein time we have rE = c2M22Mrei£E, which leads 
to the following mass estimate for the lens: 

M = 0.U M (j^) =0.14t° ;£° M (10) 

where we neglect all sources of error except the uncertainty in the Einstein 
time fit. 



6 Discussion and Summary 

We have presented the results of a search for erupting objects in the field of 
the globular cluster M22. A new cataclysmic variable, CV2, which underwent 
a DN type superoutburst in 2000 was detected. Prominent superhumps with a 
period of 128 minutes were observed in the light curve during three nights. CV2 
has an X-ray counterpart detected by the XMM-Newton telescope. The cluster 
membership of the object remains an open issue. Recent dynamical models 
of GCs (Heggie and Hut 2003) predict that close binaries should sink to the 
centers of the clusters due to encounters with neighboring stars. CV2 is located 
at a distance of 3.9 core radii from the center of M22. However, we note that 
there are several examples of close binaries located in outskirts of their parent 
clusters (eg. Bassa et al. 2003; Margon et al. 1981). The position of CV2 on the 
color-magnitude diagram of M22 and its observed X-ray luminosity are both 
consistent with the hypothesis that the variable is a cluster member. This issue 
can be easily settled by measuring the systemic radial velocity of the variable. 
The cluster itself has V ra d = — 149 km/s which differentiates its stars from the 
background population. We also detected two DN outbursts for the previously 
identified variable CV1. 

The search for erupting objects also led to the detection of a probable mi- 
crolensing event in M22. It is located 2.3 arcmin from the cluster center and 
had an amplitude of 0.75 mag. We fitted the light curve to a 5-parameter model 
of a single lens event. The most likely geometry of the event places the source 
in the Galactic bulge and the lens in the cluster. 



12 



Acknowledgements. 

PP would like to thank Arck Olcch for helpful comments during data anal- 
ysis. This work was supported by the following grants from the Ministry of 
Science and Information Society Technologies, Poland: 1 P03D 024 26 (to PP), 
1 P03D 001 28 (to JK) and 2 P03D 016 24 (to MJ). 

REFERENCES 

Alard, C, and Lupton, R.H. 1998, Astrophys. J., 503, 325. 
Albrow, M.D., de Marchi, G., and Sahu, K.C. 2002, Astrophys. J., 579, 660. 
Anderson, J., Cool, A.M., and King I.R. 2003, Astrophys. J., 597, L137. 
Baskill, D.S., Wheatley, P.J., and Osborne, J. P. 2005, MNRAS, 357, 626. 
Bassa, C.G., Vcrbunt, F., van Kerkwijk, M.H., and Homer, L. 2003, Astron. Astrophys., 409, 
L31. 

Bond, LA., et al. 2005, Astrophys. J., 620, L103. 
Dolphin, A.E. 2000a, P.A.S.P., 112, 1383. 
Dolphin, A.E. 2000b, P.A.S.P., 112, 1397. 

Fabian, A.C., Pringle, J.E., and Rees, M.J. 1975, MNRAS, 172, 15. 
Gaudi, B.S. 2002, Astrophys. J., 566, 452. 

Gendre, B., Barret D., and Webb, N. 2003, Astron. Astrophys., 403, Lll. 
Han, Ch. and Gould, A. 1996, Astrophys. J., 467, 540. 

Hannikaincn, D.C., Charles, P.A., van Zyl, L., Kong, A.K.H., Homer, L., Hakala, P., Naylor, 

T., and Davies, M.B. 2005, MNRAS, 357, 325. 
Harris, W.E. 1996, Astron. J., 112, 1487. 

Heggie, D., and Hut, P. 2003, The Gravitational Million-Body Problem - a Multidisciplinary 
Approach to Star Cluster Dynamics, Cambridge Univ. Press, Cambridge. 

Hcinke, CO., Grindlay, J.E., Edmonds, P.D., Cohn H.N., Lugger, P.M., Camilo, F., Bogdanov, 
S., and Freire, P.C. 2005, astro-ph/0503133. 

Hut, P., and Paczyhski, B. 1984, Astrophys. J., 284, 675. 

Hut, P., et al. 1992, P.A.S.P., 104, 981. 

Hut, P., et al. 2003, New Astron., 8, 337. 

Kaluzny, J., Pictrukowicz, P., Thompson, LB., Krzeminski, W., Schwarzenbcrg-Czcrny, A., 

Pych, W., and Stachowski, G. 2005, MNRAS, 359, 677. 
Landolt, A. U. 1992, Astron. J., 104, 340. 

Margon, B., Downes, R.A., and Gunn, J.E. 1981, Astrophys. J., 247, L89. 

Monaco, L., Pancino, E., Ferraro, F.R., and Bellazzini, M. 2004, MNRAS, 349, 1278. 

Paczyhski, B. 1986, Astrophys. J., 304, 1. 

Peterson, R.C., and Cudworth K.M. 1994, Astrophys. J., 420, 612. 
Richter, P., Hilker, M., and Richtlcr, T. 1999, Astron. Astrophys., 350, 476. 
Ritter, H., and Kolb, U. 2003, Astron. Astrophys., 404, 301. 

Sahu, K.C, Cascrtano, S., Livio, M., Gilliland R.L., Panagia, N., Albrow, M.D., and Potter, 

M. 2001, Nature, 411, 1022. 
Schlegel, D.J., Finkbeiner, D.P., and Davis, M. 1998, Astrophys. J., 500, 525. 
Schwarzenbcrg-Czerny, A. 1996, Astrophys. J., 460, L107. 
Skillman, D.R., and Patterson, J. 1993, Astrophys. J., 417, 298. 
Stetson, P.B. 1987, P.A.S.P., 99, 191. 

Stolz, B., and Schoembs, R. 1984, Astron. Astrophys., 132, 187. 
Trager, S.C, King, I.R., and Djorgovski, S. 1995, Astron. J., 109, 218. 
Warner, B. 1976, Observatory, 96, 49. 

Warner, B. 1995, Cataclysmic Variable Stars, Cambridge, Cambridge University Press. 
Webb, N.A., Serre, D., Gendre, B., Barret, D., Lasota, J. -P., and Rizzi, L. 2004, Astron. 

Astrophys., 424, 133. 
Wozniak, P., and Paczyhski, B. 1997, Astrophys. J., 487, 55. 
Wozniak, P.R. 2000, Acta Astron., 50, 421. 



